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ABSTRACT 

We present 7 sequential weekly observations of NGC 5548 conducted in 2007 
with the Suzaku X-ray Imaging Spectrometer (XIS) in the 0.2-12 keV band and 
Hard X-ray Detector (HXD) in 10-600 keV band. The iron Ka line is well 
detected in all seven observations and K(3 line is also detected in four observations. 
In this paper, we investigate the origin of the Fe K lines using both the width of 
the line and the reverberation mapping method. 

With the co-added XIS and HXD spectra, we identify Fe Ka and K/3 line at 
6.396lg;oo? keV and 7 - 08 -ao5 keV > respectively. The width of line obtained from 
the co-added spectra is 38^ eV (FWHM = 4200±2ooo km / s ) which corresponds 
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to a radius of 20^° light days, for the virial production of 1.220 x 10 7 M in 
NCG 5548. 

To quantitatively investigate the origin of the narrow Fe line by the rever- 
beration mapping method, we compare the observed light curves of Fe Ka line 
with the predicted ones, which are obtained by convolving the continuum light 
curve with the transfer functions in a thin shell and an inclined disk. The best- 
fit result is given by the disk case with % = 30° which is better than a fit to a 
constant flux of the Fe K line at the 92.7% level (F-test). However, the results 
with other geometries are also acceptable (P>50%). We find that the emitting 
radius obtained from the light curve is 25-37 light days, which is consistent with 
the radius derived from the Fe K line width. Combining the results of the line 
width and variation, the most likely site for the origin of the narrow iron lines 
is 20-40 light days away from the central engine, though other possibilities are 
not completely ruled out. This radius is larger than the H/3 emitting parts of 
the broad line region at 6-10 light days (obtained by the simultaneous optical 
observation), and smaller than the inner radius of the hot dust in NGC 5548 (at 
about 50 light days). 

Subject headings: galaxies : active - galaxies : individual (NGC 5548) - galaxies 
: Seyfert - quasars: emission lines - X-rays: galaxies 

1. Introduction 

The neutral, photoionized iron K emission line commonly found in the X-ray spectra of 
active galactic nuclei (AGNs) can be decomposed into a narrow core around 6.4 keV and a 
broad redshifted component, though sometimes only one of them is present (Fabian et al. 
2000; Yaqoob et al. 2001; Nandra et al. 1997, 2007). The broad, relativistic component 
could be modelled as being emitted from the vicinity of the central black hole and used 
to determine the parameters of the black hole (e.g. Brenneman & Reynolds 2006). An 
alternative explanation of this feature is the complex absorption origin (Miller, Turner & 
Reeves 2008). However, the origin of the narrow component is still unclear, but could be 
critical to understanding the central structure of AGNs. The narrow width (FWHM ~ several 
thousand km/s or less, e.g. Yaqoob & Padmanabhan 2004) suggests an origin well outside 
the continuum producing region. One of the proposed sites is the putative 'obscuring torus' 
(e.g. Nandra 2006) at ~0.1-1 pc. The evaporation radius of the dust can be estimated by 
r = 1.3L\(y 46^~i5oo 8 P c (Barvainis 1987), where L U v,46 ^ s the ultraviolet luminosity in units of 
10 46 ergs/s, and T 150 o is the evaporation temperature in units of 1500 K. For NGC 5548, since 
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Luv,46 ~ 0.01 and assuming T 150 o = 1, the evaporation radius is about 0.1 pc. However, the 
broad emission lines have similar line widths, so the more compact broad line region (BLR) 
is another possible location which cannot be ruled out simply. A BLR origin is supported by 
the quasi-simultaneous optical spectroscopic observation with Chandra observation of NGC 
7213 (Bianchi et al. 2008), which shows consistent Fe K and H/3 line widths, and by the 
rapid N H changes seen in several AGNs (Elvis et al. 2004; Puccetti et al. 2007; Risaliti et 
al. 2002), which require a BLR-like radius for the N H > 10 23 cm~ 2 absorbers. Seen from 
another angle, these absorbers must re-emit in Fe K. 

The lag between the variation of the flux of the continuum and the line can be used 
to measure the location of an emission line region (Blandford & McKee 1982) and this 
'reverberation mapping' methodology has been applied to the optical and UV broad emission 
lines (BELs) with great success (Peterson et al. 2004). However, until now the 10% or greater 
error on the flux of the 6.4 keV Fe K line (compared with the usual l%-5% error on the 
flux of BELs), and the low sampling frequency of the X-ray observations, have made it 
hard to apply this method to determine lags for Fe K lines, especially on relatively short 
timescale expected (~10 days for a BLR origin). For example, Chiang et al. (2000) found 
the flux of the iron K line in NGC 5548 was consistent with being constant using the four 
simultaneous observations of ASCA and RXTE within 25 days and another observation 
after about half a year. However, some response of the Fe K line to the continuum changes 
has been found. Markowitz et al. (2003) tried this method, analyzing the long-term RXTE 
spectra for seven sources. Although they found no evidence for correlated variability between 
the line and continuum, comparable systematic long-term (~3-4 years) decreases in the line 
and continuum were present in NGC 5548. 

We report here on a series of seven sequential X-ray observations of NGC 5548 by 
Suzaku, spaced roughly weekly. This observing campaign was designed, among other goals, 
to constrain the flux of the iron line to about 10% in each observation. This allows us, for 
the first time, to apply the reverberation mapping technique to Fe K to try to distinguish 
different geometries of the Fe K emitting region. NGC 5548 is the source best studied by the 
optical reverberation mapping technique (e.g. Peterson & Wandel 1999) and so has the best 
determined radial structure. Therefore, it is easy to determine the relative location of the 
Fe K line emitting region. In previous observations, only the Fe Ka line was detected. As 
we will present in this paper, the Kj3 line is also well detected in the Suzaku observations. 
Although very weak in some observations, Kf3 is useful to constrain the ionization state of 
iron. 

The iron line in NGG 5548 has been observed by several X-ray satellites: ASCA spectra 
suggested a relativistic broad iron Ka line in NGC 5548 with a = 3401^0 e V (68% error 



-4- 



for 4 interesting parameters, Nandra et al. 1997) and o ~ 400 eV (Chiang et al. 2000). 
However, less than two years later, only the narrow Ka line was detected in a Chandra 
observation with much higher energy resolution (38 eV vs 160 eV) but lower S /N (Yaqoob et 
al. 2001). An XMM- Newton EPIC CCD observation confirmed the Chandra result (Pounds 
et al. 2003). 

In §2, we describe the observations and the procedure of the data reduction. In this 
paper we intend to investigate the origin of the Fe K line using the reverberation mapping 
method and the width of the line. Therefore, in §3, we first fit the spectra of each observation 
to determine the flux of the continuum and the iron line. We then fit the co-added spectra 
in §3.2 to determine the mean parameters, especially the width of the iron line. In §4.2 and 
§4.3, we calculate the transfer functions in different geometries and investigate the possible 
emitting region of iron line. In §5.1, we discuss the possible origin of the iron line. In §5.3, 
we briefly discuss the implications of the intensity and the equivalent width of iron line. 
In §6 we give our conclusions. Throughout this paper we adopt the redshift of NGC 5548 
obtained from 21 cm H I measurements, i.e. z=0. 017175 (de Vaucouleurs et al. 1991). The 
errors quoted in this paper correspond to 90% confidence level (A% 2 = 2.706, Avni 1976) if 
not otherwise specified. 



2. Observations and data reduction 
2.1. Observations 

During 2007 June to 2007 August, NGC 5548 was observed by the CCD X-ray Imaging 
Spectrometers (XIS 0, 1, and 3) in 0.2-12 keV band (Koyama et al. 2007) and by the Hard 
X-ray Detector (HXD) in 10-600 keV band (Takahashi et al. 2007) on Suzaku (Mitsuda et 
al. 2007) seven times for 28.9 ks-38.7 ks each. We denote these as observations 1-7. The 
details of the observations are summarized in Table 1. 



2.2. Data reduction 

Following the standard procedures outlined in the "Suzaku Data Reduction (ABC) 
Guide (version 2)0", we used the updated Charge Transfer Inefficiency (CTI) calibration 
(Suzaku XIS CALDB 20081110) and screened the events using the xispi (Ftools 6.5) and 



1 http://heas arc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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xselect scripts provided by Suzaku team (we adopted the standard criterion in xis.event.sel 
and xis_mkf.sel)c|, respectively. 

X-ray spectra were extracted using xselect§| from all the XISs with a circular ex- 
traction region of radius 260 arcsec centered on NGC 5548 (a =14h 17m 59.5s, 5 =+25d 
08m 12.4s, J2000). Background spectra were obtained from a larger annulus around the 
source (but avoiding the calibration sources on the corners of the chips). Response matrices 
(rmf ) and effective area (arf ) files were generated with the xisrmf gerjj and xissimarf ger@ 
(estepf ile=dense and mim_photon=300000), respectively. We then combined the spectra, 
background, rmf, and arf files from the two front-side illuminated CCDs, XIS and 3, for 
each observation with addasc aspect (we denote the combined spectra by 'XIS03'). The 
spectra of XIS 1 are considered separately, as this detector uses a back-side illuminated 
CCDs. 

NGC 5548 was detected by the silicon diode PIN instrument of the HXD, but below the 
sensitivity limit of the GSO crystal scintillator instrument. We downloaded the tuned non- 
X-ray background (NXB, version 2.0) files from the Suzaku Guest Observer Facility (GOF)Q 
and then merged the good-time interval (GTI) of the NXB with that of the screened event 
files to produce a common GTI using mgtime. Then the spectra of the source and the NXB 
were extracted by xselect using the common GTI. And the dead time of the source spectra 
was corrected by hxddtcor. Since the event rate in the PIN background event file is 10 
times higher than the real background to suppress the Poisson errors, the exposure time of 
the spectra of NXB was increased by a factor of 10. The cosmic X-ray background (CXB) 
was not taken into account in the NXB file, therefore we also added a CXB component in 
the spectral fitting using the model given by the ABC guide (Section 7.3.3). The response 
file ae_hxd_pinxinome3_20080129 . rsp was used for observation 1-5, while the response file 
ae_hxd_pinxinome4_20080129 .rsp was used for observation 6 and 7 due to the changes in 
instrumental settings of Suzaku during different observation epochs. The PIN spectra were 
multiplied by a constant cross-normalization factor of 1.16 to account for the differences in 



-http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/xisrepro.xco, http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/xis_e 



http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/xis_mkf.sel 



3 http: / /heasarc. nasa.gov/docs / software /lheasoft /ftools /xselect/index.html 
4 http: / /heasarc. nasa.gov/docs/suzaku/analysis/xisrmfgen. html 
5 http://heas arc. gsfc.nasa.gov/docs/suzaku/analysis/xissimarfgen/ 
6 http:/ /heasarc. gsfc.nasa.gov/lheasoft/ftools/fhelp/addascaspec.txt 
7 ftp:/ /legacy, gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.0_tuned/ 
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calibration between XIS and PINH. 



3. Spectral fitting 

To perform the reverberation mapping calculation, we should first determine the flux of 
the continuum and the Fe K line. In order to avoid the influence of the complex absorption 
below 3 keV (Detmers et al. 2008; Steenbrugge et al. 2003), we only analyze the XIS spectra 
in 3-10 keV band. The warm absorber seen in the <3 keV spectra is presented by Krongold 
et al. (2009). In this paper, we will focus on the property and origin of the iron emission 
line. A global fit to the entire Suzaku spectral band for the entire campaign will be presented 
in a subsequent paper. 

For each observation, we fitted the spectra of XIS03 simultaneously with the spectra of 
XIS 1 in XSPEC (version 12.4). Although both of the energy resolution and the effective 
area of XIS 1 are lower than XIS 0, 3 in 4-10 keV band (the background level of XIS 1 is 
also higher than XIS 0, 3 in this band), it is still useful to reduce the error of the intensity of 
the Fe K line. We fixed the Galactic absorption column at 1.63 xlO 20 cm -2 (Murphy et al. 
1996), and fitted the spectra with a single power law. The continua were well described by 
a single power law except for the region of the iron Ka and K/3 lines around 6.4 keV and 7.0 
keV, respectively (see Figure 1). Therefore, we added two Gaussian lines to describe these 
features. It is important to construct a self-consistent model to describe all the components 
of the continuum and theoretically predict the strength of the Fe K line, as in Murphy & 
Yaqoob (2009). However, this is not the purpose of this paper. Our purpose is only to 
simply model the Fe K line using a Gaussian line and determine the flux and width of it 
and then to perform the reverberation calculation. A Ka line is required by all the seven 
spectra (>6 a), while K/3 line is only required by four observations (>2 a for 1, 2, 4, and 
7). The 90% upper limit of the flux of the K/3 line is determined in observation 3, 5, and 
6. Since the K/3 line is weak, we fixed the width of K/3 to be the same as that of Ka. Fe 
K/3 line was also detected in other sources, e.g. NGC 2992 (Yaqoob et al. 2007) and Mrk 
3 (Awaki et al. 2008). The detailed result of the fitting is given in Table 3, though we are 
mainly concerned about the flux of the continuum and the intensity of the Fe K line. Due to 
the weak reflection component in NCG 5548 (e.g. Pounds et al. 2003), if we simultaneously 
fit the continua in the XIS and PIN spectra of each observation using the pexrav model 
instead of the power law model, the intensity of the Fe K line will systematically decrease by 
10% for each observation (the strength of the reflection component R ~ 0.2 — 1.5). However, 



http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/watchout.html 
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this change will not influence the reverberation mapping result in §4.3, since it degenerates 
with the normalization of the transfer function (see the detail of the transfer function in 
§4.2). The PIN data is also not helpful to reduce the error of the flux of the Fe K line 
due to the additional uncertainty introduced by the reflection component and large error in 
the PIN spectra. The detailed discussion about the variation of the X-ray continuum and 
simultaneous UV/optical data will be presented in another paper. 

3.1. Continuum light curve 

Having fitted the spectra, we calculated the observed flux of the continuum in the 3-10 
keV band. The result light curve is shown in Figure 3a. To better sample the continuum 
light curve, we also utilized one observation of the continuum from the simultaneous Swift 
campaign. Since the observation times and the flux of other Swift data are quite similar to 
that of the Suzaku data, we will not include them in this paper. The details of the Swift 
campaign will be discussed in Grupe et al. (2009). 28 observations of the continuum with 
PCA on RXTE before the Suzaku campaign are also used, since we are looking for the lag 
between the variation of continuum and line. However, due to the short exposure time, we 
cannot determine the flux of the line in any of these additional observations. The details of 
the RXTE and Swift observations are summarized in Table 2. 

3.2. Narrow iron lines 

To determine the mean parameters of the Fe Ka and K/3 lines more accurately, especially 
the width of the line, we added the spectra from the seven observations together using 
addascaspec (the XIS03, XIS1, and PIN spectra were added separately and then fitted 
simultaneously). The net counts of the source in 3-10 keV band in the total XIS03 and XIS1 
spectra are 160349 and 75332, respectively. The net counts of the source in 12-35 keV band 
in the total PIN spectra is 22243. We found that the XIS1 spectra could not provide any 
useful constraint on the width of the Fe K line, since the a of the Gaussian line is pegged 
at eV (the 90% upper limit is 30 eV). Therefore, we will only utilize the co-added XIS 03 
and PIN spectra to determine the width of the Fe K line. 

If we only fitted the co-add XIS03 spectra using the model in §3, i.e. a power law and 
two Gaussian lines, the width is a = 50tJs eV. However, since the weak reflection component 
could influence the width of the Fe K line, we then simultaneously fitted the co-added XIS 03 
and PIN spectra using the pexrav model and two Gaussian lines. The derived parameters 
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are given in Table 4, where the strength of the reflection component (R ~ 0.8) could explain 
the Fe K line. As shown in Figure 1, only a narrow iron line is clearly present in the spectra, 
with a width, a =38^}g eV. This value is consistent with previous results: i.e. a =41^4 eV 
obtained by HEG+MEG on Chandra (Yaqoob et al. 2001) and a =40^ eV (MOS) and 
64±24 e V (PN) obtained by XMM-Newton (Pounds et al. 2003). 

Since the peak energies of the line and the parameters of the continua are somewhat 
different for each observation, we simultaneously fitted the spectra of all observations in 
XSPEC to test whether the line was broadened artificially in the co-added process. We 
required the width of line to be the same for all observations and kept other parameters free. 
The obtained width is only slightly smaller than that from the co-add spectra by about 2 
eV, which implies the co-added method has not significantly broadened the width. 

The width obtained by the co-add spectra is inconsistent with zero at 2.2 a and corre- 
sponds to FWHM= 4200 km/s and a radius of 5.2xl0 16 cm or 5.2xl0 3 r g (r g = GM/c 2 ), for 
the 6.71xlO 7 M black hole in NGC 5548 (Peterson et al. 2004). In §5.1, we will estimate 
the location of the emitting material of the line and discuss the origin of the line. 

To test the presence of the broad, relativistic Fe Ka line found by ASCA, we tried the 
diskline model in XSPEC to fit the Ka line. We fixed the inner and outer disk radii at Qr g 
and 1000r 9 (r- g = GM/c 2 ), respectively. The index of the power law emissivity was frozen 
at -2.5 (the mean value obtained in Nandra et al. 1997). If we thaw the index, it will be 
pegged at the positive upper limit in XSPEC, i.e., the flux of the line is dominated by the 
outer disk and therefore it is not a disk line at all. It was found that the diskline model 
cannot describe the Ka line alone, since x 2 was higher by 128 than for a narrow line with 
the same number of parameters. The same conclusion was also obtained by Yaqoob et al. 
(2001). 

Next, we investigated the result if we fit the Ka line with a disk line and a Gaussian 
line. Besides the constraint on the diskline model mentioned above, we also fixed the value 
of the inclination angle at 31 degrees, which is the best-constrained value from the ASCA 
data (Yaqoob et al. 2001). We found the central value of the intensity of the disk line 
is pegged at and the 90% upper limit is 4 x 10~ 6 photons/cm 2 /s. Therefore, any broad 
component must be > 5 times weaker than the narrow component and we will not include 
this component in the following discussion. 

We show the confidence regions for the peak energies of Ka and K/3 in Figure 2a. The 
expected values of different iron ionization states (Palmeri et al. 2003; Mendoza et al. 2004; 
Yaqoob et al. 2007) are also shown. Since the peak energy may be influenced by the residuals 
in the energy scale calibration, we extracted the spectra of the 55 Fe calibration sources on 
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the corners of the CCDs and added them together (see the inset in Figure lb). Using a 
Gaussian line to fit the Mn Ka line, we found the peak energies are 5.8921q;qq} keV and 
5.893^q;oo^ keV for XIS03 and XIS1, respectively. The expected value of the Mn Ka line is 
5.895 keV. This result is well within the accuracy of the absolute XIS energy scale given by 
the Suzaku Technical Description as O.2%0. The fit widths of Mn Ka line are 12jJ eV 
and 16lg eV for XIS03 and XIS1, respectively. Therefore, the peak energy obtained by the 
spectral fitting is reliable and we could conclude that the iron emission line in NGC 5548 is 
most likely to be dominated by relatively low ionization states, <Fe XIII (99% confidence). 
Actually, the observed Fe K line could be a blend of several ionization states. 

The confidence regions of the intensities of the Fe Ka and K/3 lines are shown in Figure 
2b. The expected ratio of I(K/3) to I(Kcn) varies from 0.12 to 0.20 for different ionization 
states (dashed lines in Figure 2b, Palmeri et al. 2003; Mendoza et al. 2004;). Our result is 
fully consistent with the expected value. However, due to the weakness of K(3, the error is 
still too large to constrain the ionization state precisely. 

4. Time variability analysis 

4.1. Fe Ka line light curves 

The light curve of Fe Ka line is shown in Figure 3b each being determined to ~10%. 
As K/3 is weak in some observations, we will only discuss the result for Ka below. As shown 
in Figure 3a, the flux of the continuum changed strongly during the seven observations. The 
highest value (2.65 xlO -11 ergs/cm 2 /s, for observation 5) is about 4 times that of the lowest 
(6.89xl0~ 12 ergs/cm 2 /s, for observation 1). However, perhaps due to the relatively large 
error, the flux of the line is consistent with being constant, Xmin =3.05 (P=80%) and the 
value of the corresponding line constant flux is 2.24xl0 -5 photons/cm 2 /s. 

A constant Fe Ka flux is the simplest solution to the Fe Ka light curve, and it requires 
the emitting region should be far from the central engine to smooth out the variation in 
short time scale. According to the result in §4.3 and Figure 6b, it should be larger than 100 
light days. However, as shown in Figure 1 and 2 in Markowitz et al. (2003), the presence 
of the variation of the Fe K flux in time scale shorter than 100 days indicates the emitting 



9 http:/ /www. astro. isas.jaxa.jp/suzaku/doc/suzaku_td/ 

10 We found that using the calibration taken for the same time interval was important. A preliminary 
analysis using the calibration available when the observations were made showed a 20 eV offset, which was 
puzzling. The amplitude of the Mn Ka apparent energy variation was also about 20 eV 
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region must be smaller than 100 light days. The emitting region required by the constant flux 
fitting is also inconsistent with that required by the line width (see §5.1). Therefore, we have 
investigated whether there is a better and more self-consistent solution than the constant 
flux fitting. To quantitatively access the emitting location of the Fe Ka line, we should 
perform the reverberation calculation. The cross-correlation function (CCF) is usually used 
in the optical reverberation mapping to determine the lag between the continuum and the 
emission line. However, due to the very few data points of the light curve of the Fe K line 
and the large error of them, we did not find any significant peak in CCF. Therefore, we will 
approach this problem in a different way. To do so we first calculate the transfer function 
for Fe Ka in 84.2. 



4.2. Transfer function 

As in the case of the optical broad emission lines, we calculate a transfer function for 
Fe Ka, i.e. the response of the flux of the line to a 5 function change in the continuum. We 
consider two cases: (1) a spherical region and (2) an inclined disk. 

(1) We assume the emitting material is spherically distributed around the center with 
an inner (r min ) and outer (r max ) radius. The transfer function in this model is then simply a 
constant between and 2r min /c, and then decays to at 2r max /c (see Figure 4a and Peterson 
1993). 

(2) For a thin disk, the general form of the transfer function is two-peaked (Welsh & 
Home 1991). Figure 4b shows the result for different values of the inclination angle % (the 
angle between the normal of the disk and the line of sight). 

The shape of the transfer function depends on the form of the "responsivity" e(r) (see 
Figure 4a), which combines the effects of the distributions of the number density of clouds 
and the emissivity per cloud. We assumed that the responsivity is a power law e(r) cc r~ a 
and the normalization is adjusted to fit the observed light curve. The power law form is 
simple and somewhat arbitrary. However, since we only consider the thin shell case (i.e. 
Ar <C r, or equivalent to the locally optimized clouds model), the result is not sensitive to 
the detailed form of the responsivity nor the index of the power law (see Figure 6a). a = 3 
is adopted in the calculation in §4.3. 

We could then obtain the predicted light curve of the line by convolving the transfer 
function with the light curve of the continuum flux. In the thin spherical shell and disk 
cases discussed in §4.3, it can be proved that the the lag time, r, obtained by the CCF just 
corresponds to the radius of the emitting region, i.e. r = r/c (Koratkar & Gaskell 1991). 
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4.3. Comparison with the observed light curve 

Since the observation data of the continuum are still too few to produce a complete light 
curve, we performed linear interpolation between data points to convolve the light curve with 
the transfer function. 

We used the light curve of the continuum in 3-10 keV band in order to have precision 
measurements. Similar bands (e.g. 2-10 keV) have been adopted in the previous attempts to 
determine the lag between the variation of the continuum and the Fe K line (e.g. Markowitz 
et al. 2003, 2009), although only photons from above the ionization threshold can actually 
lead the emission of an Fe K photorJ"1. As shown in Figure 5, the flux in 3-6 keV band is 
tightly and nearly proportionally correlated with the flux in 8-10 keV band. Any possible 
effects due to rapidly changing Nh (Risaliti et al. 2005) must therefore be small. As a 
result, the results of the following calculations will not be sensitive to the adopted energy 
band of the continuum. The change of the continuum slope is also related to the total 
amount of the ionization flux. However, this is a minor factor compared with the change 
of the normalization of the continuum. Specially, this effect should be even small for NGC 
5548, since e.g. Sobolewska and Papadakis (2009) showed that the continuum slope is nearly 
independent on the flux. We find the same lack of dependence (Table 3). 

We fixed the width of the emitting region at 0.1 light days, substantially smaller than the 
likely radius of the Fe K emitting region, and varied r m j n . After convolving the light curve of 
the continuum with the transfer functions in different geometries and r min , we compared the 
predicted light curve of the line with the observed one to find the minimum value of x 2 (since 
the error on the flux is asymmetric, we conservatively adopt the larger one to calculate x 2 )- 
Figure 6b shows x 2 vs r mm- All the curves show pronounced minima (somewhat surprisingly, 
given the weak structure in the Fe Ka light curve) and finally decrease towards the result of 
the constant fitting (the horizontal dashed in Figure 6b) with increasing r min , as the result of 
the significant smoothing effect with large r min . The values of minimum x 2 m the spherical 
thin shell case and disk case with i = 15° and i = 30° are well below that of the constant 
fitting, and we will only discuss these three cases below. With the smallest value of x 2 , the 
best-fitting light curve in the disk with i = 30° is improved at 92.7% level compared with 
the constant fitting. We show the predicted light curves corresponding to the minima of x 2 
in these three cases in Figure 7. Due to the few data points, we cannot tightly constrain the 
inclination angle. Except for the the disk case with i = 15°, the positions of the minimum 
X 2 in the spherical case (27 l 22 light days) and the disk case with i = 30° (37 light days) 



n A similar approximation is used in reverberation mapping of AGNs (e.g. Peterson et al. 2002; Bentz et 
al. 2007) 
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are similar, and they are smaller than the inner radius of the dust (47-53 light days, see the 
discussion in §5.1). The value of the best-fitting r min in the disk case with i = 15° is 61^ 
light days, which is beyond the inner radius of the dust of NGC 5548 (such optically thick 
region will significantly absorb the photons of Fe K lines) and only marginally consistent 
with the small tail of the 90% confidence interval of the emitting region inferred from the 
widths of Fe Ka line (see the error bar in Figure 9 and it should be noted that the error bar is 
quite asymmetric). In addition, % = 15° is also well smaller than the best-fitting inclination 
angle of the disk line in the ASCA (i = 31.5° ± 6.5°, Yaqoob et al. 2001), though it is not 
very reliable, due to the absence of the broad component of the Fe K line in the following 
Chandra, XMM-Newton, and our Suzaku observations. Therefore, the disk case with i = 15° 
is quite unlikely, and we will only focus on the spherical case and the disk case with i = 30° 
when discussing the origin of the narrow Fe Ka line in detail in §5.1. 



Assuming the geometry and dynamics of the emitting region of the Fe K line are the 
same as that of the H/3 line, and using the virial relation for the H/3 line, a 2 r/G = 1.220 x 10 7 
M Q (Peterson et al. 2004), the radius of the Fe K emitting region can be derived using the Fe 
Ka width, cr, obtained in §3.2. The derived radius using the width obtained by the co-added 



In Figure 8, we plot the line width size against r m j n from the reverberation analysis 
(Figure 6b) for both the spherical thin shell case and the disk case with i = 30° with 90% 
confidence intervals for both quantities. The disk case with i = 30° and the spherical case 
are consistent with the width of the co-added Fe Ka line Combining the above results, 
we conclude that the origin of the narrow iron line in NGC 5548 is likely to be ~ 20 — 40 
light days away from the continuum source, for the geometries considered. However, the 



12 Since the virial production is an observed quantity and the radius derived from width will be compared 
with that also obtained from the reverberation mapping method, no additional geometry factor is required. 

13 Thc calibration residual in the width of Mn Ka line derived in §3.2 (12±| eV and 16±f eV for XIS03 
and XIS1, respectively) is partly due to the systematic error on the calibration of the non-Gaussian response 
function of XIS (Koyama et al. 2007). The energy resolution at the center of the CCD chip is also slightly 
better than that at the corner, but this difference is smaller than the systematic error on the calibration. 
Therefore, the true width of the Fe K line could be smaller than the observed one by a few eVs due to the 
above factors. However, these effects could not be accurately corrected simply. 



5. 



Discussion 



5.1. Location of the Fe K emitting region 
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other possible origins are not completely ruled out due to the incompleteness of the light 
curves, our model-dependent method, and the sizeable error on the width of the Fe K line. 

Since NCG 5548 is one of the best-studied AGNs with reverberation mapping, the 
locations of the BELs are well-determined, especially for the H/3 line. The lag between the 
flux of H/9 line and the continuum varies from 6.5 light days to 26.5 light day depending 
on the luminosity of the continuum (Bentz et al. 2007; Cackett & Home 2006). Since 
the correlation between the lag time and the continuum flux is more significant than the 
correlation between the lag time and the width of H/3 line (Bentz et al. 2007), and the broad 
component of the H/3 line is weak during the Suzaku campaign, we will utilize the continuum 
flux to predict the radius of the H/3 line region. 

From the simultaneous optical spectra of NGC 5548 from FLWO FAST spectrograph 
(June 19-23, 2007), we measured the monochromatic flux at 5100 A and used this flux to 
estimate the radius of the H/3 emitting region at the time of the Suzaku campaign. After 
calibration with a standard star, the optical spectra were put on an absolute calibration scale 
assuming a constant flux of F([OIII] 5007) g ^ anc ^ arc j = 5.58xl0~ 13 ergs/s/cm 2 (Peterson et 
al. 1991) and found the monochromatic flux at 5100 A, F(5100 A) ]- )Serve( j = 5.13xl0~ 15 
ergs/s/cm 2 /A. For the aperture size of 5 arcsec x 7.5 arcsec, the contribution of the host 
galaxy at 5100 A is 4.47xl0" 15 ergs/s/cm 2 /A (Bentz et al. 2006). To account for the 
effect of the aperture and the difference between telescopes, the flux measured by FLWO 
should be converted by the coefficients given in Peterson et al. (2002), i.e. F(5100 A) true = 
<^F(5100 A.) k servec [ — G. We averaged the results using the different coefficients (i.e. (p and 
G) determined in years 8 and 9-13 of the campaign (see details in Peterson et al. 2002), and 
then excluded the flux of the host galaxy (4.47xl0~ 15 ergs/s/cm 2 /A, Bentz et al. [2006]) . 
The final 5100 A flux of the AGN in NGC 5548 is (1.7 ± 0.7) x 10~ 15 ergs/s/cm 2 /A during 
the Suzaku campaign. With the relation between the emitting region of H/3 line and F(5100 
A) (Figure 5 in Bentz et al. [2007]), we found r(H/3)= 8 ± 2 light days during the Suzaku 
observations. 

From the result of dust reverberation of NGC 5548, we also know the inner radius of 
the hot dust in this source is 47-53 light days (from the lag of K-band relative to V-band, 
Suganuma et al. 2006). 

We plot the radial locations of Fe Ka, R/3 and the inner dust in Figure 9 (note that this 
is a linear plot). It is likely that the location of the Fe K emitting region is closer in than the 
inner radius of the dust, but slightly farther out than H/3, i.e. the origin of the narrow iron 
line lies in the outer part of the BLR [(0.5 — 1.0) x 10 4 r g \. This result is consistent with the 
region emitting the optical Fe II emission lines, which is less than several light weeks from the 
continuum emitter (Vestergaard & Peterson 2005). This region is also consistent with the 
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intermediate line region proposed by Zhu et al. (2009) and the radius of the low ionization, 
low velocity component of the warm absorber in NGC 5548 (r<3 pc, Krongold et al. 2009). 
Despite the uncertainty of the column density of the Fe K emitting region, it is surely much 
higher (see §5.3) than that of the warm absorber found in NGC 5548 (10 21 ~ 10 22 cm -2 , 
Krongold et al. 2009). Therefore, It clearly does not lie along our line of sight. 

5.2. Physical condition 

We can use the location derived above to estimate the density of the emitting region 
from the ionization parameter and the distance. Kallman et al. (2004) investigated the 
dependence of the line profile of Ka and Kj3 on the ionization parameter. The peak energies 
of Ka and K/3 in NGC 5548 (see Figure 2a) prefer a relatively low ionization parameter, 
i.e. log£ = — 1 ± 1, where £ = L/nR 2 (erg s _1 cm), L is the Fe K ionizing luminosity of 
the continuum source (erg s -1 ), n is the density of the emitting region (cm -3 ), and R is 
the distance to the continuum source (cm). With logL = 43.7 ± 0.3, log£ = — 1 ± 1, and 
R = 30 ± 10 light days, we obtain logn = 10.9 ± 1.1, which is comparable with the density 
of the gas in the BLR of NGC 5548 (Ferland et al. 1992; Goad & Koratkar 1998; Kaspi & 
Netzer 1999). 

5.3. Theoretical intensity of the Fe line 

The intensity and the equivalent width of the iron line can be estimated theoretically 
(Krolik & Kallman 1987; Yaqoob et al. 2001). We found the column density N H > 10 23 cm" 2 
is required to produce the observed intensity and the equivalent width of the Fe Ka line in 
our Suzaku observations. However, as pointed out by Miller et al. (2009) and Yaqoob et al. 
(2009), due to the self-absorption effect and the Compton scattering, for Nh > 10 23 cm -2 , 
the relation between the intensity of the Fe K line and Nh or the abundance of iron is quite 
non-linear, and the intensity of the Fe K line also significantly depends on the geometry of 
the emitting region and the observing angle. Since the detailed investigation of the geometry, 
N H) and the abundance of the emitting region of Fe K lines is much beyond the scope of 
this paper, we will not further discuss the constraints obtained from the intensity and the 
equivalent width of the Fe K lines. 
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6. Conclusions 

We analyzed the iron Ka and K/3 lines in spectra of NGC 5548 obtained by Suzaku XIS 
and summarize our results as follows. 

1. The iron Ka line was well detected (> 6<x) in all seven observations and the Kf3 line 
was also detected (> 2a) in four observations (1, 2, 4, and 7). 

2. Only a narrow iron line was found in the spectra. The line width obtained by the 
added spectra is 38^}g eV, which is consistent with the results of Chandra and XMM-Newton. 
Assuming the same virial relation as that of the H/3 line, this width corresponds to a radius 
of 20^° light days. Any relativistically broadened disk line must be a factor of 5 weaker 
than the narrow component in flux at 90% confidence level. 

3. We compared the observed peak energies and intensity ratios of Ka and K/3 lines 
with the expected value and found they are consistent with the low ionization states of iron, 
i.e. lower than Fe XIII, at the 99% confidence level. 

4. The Fe Ka line is consistent with being constant over the 50 days of the Suzaku 
campaign, although the 3-10 keV continuum varies by a factor of 4. It is shown that a 
location at >100 light days is consistent with the data (see the discussion in §4.3 and Figure 
6b), but this is not a unique result. 

5. To further access the location of the iron lines using the light curve, we calculated the 
transfer functions in spherical and disk geometries, and compared the predicted light curves 
with the observed one. The value of x 2 is smallest in the disk case with i = 30°, which 
is better than the constant fitting at the 92.7% level. The spherical thin shell case is also 
acceptable (P=81%). The inferred emitting radii are 27^ 2 light days in the spherical case 
and 37^5 light days in the disk case with i = 30°, which are consistent with that obtained 
from the width of iron lines. 

6. Combining the Suzaku constraints, the most likely origin of the narrow iron lines is 
about 20-40 light days away from the central engine, i.e. the outer part of BLR (5.2 xlO 3 - 
1.0 xlO 4 r g ). However, we could not completely rule out other possible origins. 

The approaches used in this paper offer a valuable tool for determining the size and 
structure of the inner regions of AGNs, although we stress again this method is model 
dependent. The constraint on the emitting region of the narrow Fe K line obtained by the 
width will be greatly improved by upcoming calorimeters, which have > 10 time better the 
energy resolution, of only a few eV. If future X-ray satellites (e.g. Astro-H, IXO, and Gen-X) 
with larger effective area could reduce the error of the flux of the Fe K emission line by even 
a factor of 2, then it will be possible to distinguish different geometries from the constant 
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flux fitting. Higher sampling frequency campaign, preferably over a longer baseline, is also 
desirable to obtain a cross-correlation function with a quality comparable to or better than 
current optical observations. 
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of Education at Tsinghua University, Directional Research Project of the Chinese Academy of 
Sciences under project no. KJCX2-YW-T03 and by the National Natural Science Foundation 
of China under grant nos. 10821061, 10733010, 10725313, and by 973 Program of China 
under grant 2009CB824800. This work has been funded by CXC grant G07-8136A and 
Suzaku Grant NNX08AB81G. 
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Fig. 1. — The best fit model (solid lines) and the added spectra of XIS03 (plus). The inset 
(a) is the residual of the spectra of XIS03 to the best fit of a single power law. The inset (b) 
is the added spectra of Mn Ka line of the calibration sources on XIS03. 
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Fig. 2. — The confidence region of the peak energies (a) and intensities (b) of Ka and K/3 
obtained from the fitting of the co-added XIS03 and PIN spectra. The contours from inner 
to outer correspond to Ax 2 =1.00, 2.71, and 6.63 (68%, 90%, 99%), respectively. The crosses 
in (a) are the predicted values of the peak energies in different ionization states of iron. The 
lower and upper dashed lines in (b) correspond to the I(K/3)/I(K«) ratios of 0.12 and 0.20, 
respectively. 



- 19 - 



40 



10 



10 

aj 25 

l-M 



I • Suzaku 

I + Swift 
x RXTE 



1 

o 20 - 



It i 

1 T I 



i 

15 - 

I i 



i I t 1 I I 



10 - 



5 - 



I 



-100 -50 

Time (day) 



50 



X 2 =3-05 



10 15 20 25 30 35 40 45 50 

Time (day) 



Fig. 3. — (a) The light curve of the flux of the continuum in 3-10 keV band. The small error 
bar of the flux of the Suzaku observation is omitted, (b) The light curve of the flux of Ka 
line. The solid line in (b) is the result of the constant fitting (see the text in §4.1). 
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Fig. 4. — The transfer functions for the spherical (a) and disk cases (b). We choose r m j n = 10 
light days and r max = 20 light days for all curves, which are normalized to 1 at the maximum, 
(a) a = 2 (blue), a = 3 (red), and a = 4 (green), (b) a = 3 for all curves, i = 30° (blue), 
i = 45° (red), and i = 60° (green). 
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Fig. 5. — The top panel shows the correlation between the flux in 3-6 keV band and that in 
8-10 keV band of all seven observations. The solid line is the best-fitting straight line across 
the origin. The error bars of the fluxes are omitted, since they are smaller than the symbols. 
The bottom panel shows the ratios between the flux in 8-10 keV band and that in 3-6 keV 
band. 
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Fig. 6. — (a) Comparison with different values of a in the spherical case, a = 2 (blue), 
a = 3 (red), and a = 4 (green). Since the original curves are nearly identical, the curves 
of a = 2 and a = 4 are vertically shifted by -0.5 and 0.5, respectively, (b) The curves of 
X 2 with different r m i n for spherical and disk cases. Spherical (black), disk i = 15° (blue), 
i = 30° (red), i = 45° (green), and i = 60° (yellow). The horizontal dashed line is the x 2 °f 
the constant fitting (see the text), (c) Enlarged plot of (b). The error bar at the bottom is 
the 90% confidence interval inferred from the line width obtained by the simultaneous fitting 
of all spectra (see Figure 8). 
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Fig. 7. — The comparison between the observed and predicted light curve of the flux of the 
line with the minimum of x 2 i n different cases, i.e. (a) spherical case, (b) disk case with 
% = 15°, (c) disk case with % = 30°. 
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Fig. 8. — The relationship between a of the line and the radius of the emitting region inferred 
from the black hole mass. The vertical dashed and solid lines correspond the central value 
and the 90% confidence intervals of the line width obtained by the co-added XIS 03 and PIN 
spectra, respectively. The horizontal lines correspond to the 90% confidence intervals of the 
emitting region inferred from the light curve in the spherical case (solid) and the disk case 
with i = 30° (dashed, see Figure 6c). 
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Fig. 9. — The summary of the location of the Fe Ka line emitting region. We show intervals 
derived from the variation of the flux and the width of the Fe K line, with the locations 
of H/3 and the inner radius of the dust for comparison. The solid and dashed error bars 
correspond to 68% and 90% confidence intervals, respectively. 
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Table 1: Suzaku observation log. 



Sequence number 


Obs. ID 


Start Date & Time 


Exposure time (ks) 


3-10 keV Count Rate 
(10~ 3 photons/cm 2 /s) 


1 


702042010 


2007-06-18 UT 22:28:15 


31.1 


0.772 


2 


702042020 


2007-06-24 UT 21:53:31 


35.9 


1.29 


3 


702042040 


2007-07-08 UT 10:02:55 


30.7 


2.35 


4 


702042050 


2007-07-15 UT 13:57:39 


30.0 


1.62 


5 


702042060 


2007-07-22 UT 10:40:25 


28.9 


3.05 


6 


702042070 


2007-07-29 UT 04:20:44 


31.8 


2.04 


7 


702042080 


2007-08-05 UT 00:37:46 


38.8 


1.12 
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Table 2: RXTE and Swift observation logs. 

Obs. ID Start Date & Time 



RXTE 



92113- 


-07-40-00 


2006-12- 


■05 


UT 


03: 


:35: 


:48 


92113- 


-07-41-00 


2006-12- 


■11 


UT 


12: 


:59: 


:06 


92113- 


-07-42-00 


2006-12- 


•21 


UT 


16: 


:34: 


:27 


92113- 


-07-43-00 


2006-12- 


■23 


UT 


11: 


:08: 


Ml 


92113- 


-07-44-00 


2007-01- 


•02 


UT 


13: 


:09: 


:22 


92113- 


-07-45-00 


2007-01- 


•10 


UT 


10: 


:09: 


:27 


92113- 


-07-46-00 


2007-01- 


■14 


UT 


21: 


:15: 


:16 


92113- 


-07-47-00 


2007-01- 


•21 


UT 


16: 


:35: 


:22 


92113- 


-07-48-00 


2007-01- 


•29 


UT 


17: 


:19: 


:33 


92113- 


-07-49-00 


2007-02- 


•05 


UT 


20: 


:22: 


:13 


92113- 


-07-50-00 


2007-02- 


■13 


UT 


12: 


:12: 


:30 


92113- 


-07-51-00 


2007-02- 


•19 


UT 


08: 


:05: 


:31 


92113- 


-07-52-00 


2007-02- 


■26 


UT 


07: 


:25: 


:07 


92113- 


-07-53-00 


2007-03- 


■04 


UT 


07: 


:48: 


:42 


92113- 


-07-54-00 


2007-03- 


•13 


UT 


11: 


:30: 


:51 


92113- 


-07-55-00 


2007-03- 


■18 


UT 


18: 


:56: 


:43 


92113- 


-07-56-00 


2007-03- 


■26 


UT 


13: 


:58: 


:42 


92113- 


-07-57-00 


2007-04- 


•01 


UT 


16: 


:10: 


:48 


92113- 


-07-58-00 


2007-04- 


■09 


UT 


04. 


:12: 


:03 


92113- 


-07-59-00 


2007-04- 


•16 


UT 


03: 


:23: 


:00 


92113- 


-07-60-00 


2007-04- 


■23 


UT 


06: 


:05: 


:57 


92113- 


-07-61-00 


2007-04- 


■29 


UT 


12: 


:54: 


:16 


92113- 


-07-62-00 


2007-05- 


■07 


UT 


17: 


1 1 


:43 


92113- 


-07-63-00 


2007-05- 


■14 


UT 


11: 


:42: 


:04 


92113- 


-07-64-00 


2007-05- 


■21 


UT 


11: 


:48: 


:48 


92113- 


-07-65-00 


2007-05- 


■28 


UT 


03: 


:30: 


:13 


92113- 


-07-66-00 


2007-06- 


■04 


UT 


01: 


:59: 


:35 


92113- 


-07-67-00 


2007-06- 


■11 


UT 


07: 


:33: 


:05 



Swift 

00030022061 2007-07-01 UT20:25:01 



Exposure time (s) 3-10 keV Flux 

(10~ 12 ergs/cm 2 /s) 



919 


12.0 ± 0.9 


885 


18.1 ± 1.0 


961 


12.2 ± 0.9 


830 


15.8 ± 1.0 


914 


24.0 ± 1.0 


1259 


43.6 ± 1.0 


926 


29.3 ± 1.1 


913 


28.3 ± 1.0 


948 


21.5 ± 0.9 


1171 


20.4 ± 0.9 


895 


36.0 ± 1.1 


896 


41.9 ± 1.1 


1002 


30.2 ± 1.1 


1354 


39.6 ± 0.9 


614 


23.2 ± 1.2 


572 


26.4 ± 1.3 


994 


19.1 ± 0.9 


982 


15.7 ± 1.0 


893 


19.0 ± 1.0 


955 


14.2 ± 1.1 


961 


16.3 ± 0.9 


1041 


10.7 ± 0.8 


928 


9.4 ± 0.9 


950 


15.5 ± 0.9 


922 


13.3 ± 0.9 


931 


15.9 ± 1.0 


949 


22.7 ± 1.0 


888 


18.3 ± 1.1 



2503.18 10.5 



Table 3: The derived parameters from the spectral fits of the XIS data of each observation (with error corresponding 
to 90% confidence level except flux and intensity). 



Sequence 


Continuum 




Ka 






K/3 










number 


pa 


Flux b 


E c 


a d 


r 


EW f 


E c 


I s 


EW f 


X 2 /dof 




Av 2 ' 

^A2 


1 

2 
3 
4 


1 41+0-07 

i cc+0.05 
i - oo -0.05 
i «o+0.03 

1 co+0.04 
i - do -0.04 


6-89+-- 

11.28±8;|g 

20.25^.16 
14.2418;^ 


6.405+roi b 7 

6.384±g:8iS 

6.374lg;gS 
6-380±Hll 


461^ 
04+31 


no+0.23 

z.uz_ 24 

2 40+ - 25 
2 1fi+ - 30 

91+O.28 
z - zl -0.27 


215+^ 
155lf 6 
77±»» 

112l 4 4 8 8 


7-oo±«:gl 
7.o8i8:8^ 

7.04J 

7 05+ 04 

' - uo -0.04 


5 6+ 1 ' 7 

°-°-1.8 
re c+2.0 
0.0_ 2 

<5.5 
7 r+2.2 

'•°-2.1 


681S 
49+51 

<40 


78.9/114 
197.7/177 
355.1/338 
227.3/214 


10.2 
7.4 

11.7 


122.4 
122.8 
65.7 
98.9 


5 


1 fil+ 003 

i - oi -0.03 


26.4618;S 


6 401 + 029 

°-* ul -0.031 


fi2 +49 
DZ -62 


9 9 r +0.43 
z - zd -0.36 


6lt 3 34 


7.04J 


<5.3 


<31 


427.9/406 




54.6 


6 
7 


1 r 7 +0.03 
1 - 01 -0.03 

1 ro+0.04 

1 - oo -0.04 


17.80±g-il 

q ocr+0.20 
y - OO -0.24 


6.303±g;g« 

6.4i6l8;8g 


04+29 

90+28 
z0 -23 


9 fi 9+0.29 
z - uz -0.26 
9 1 7 +0.19 
Z - Ll -0.21 


106i 3 J 
162l| 


7.04J 
r 07+0.09 

D - y '-0.08 


<7.0 

4-o±i:g 


<50 
04+43 

— 34 


303.2/308 
205.0/216 


6.9 


124.1 
166.1 



"Photon index 

h Observcd flux in 3-10 keV band (10~ 12 ergs/cm 2 /s) with 68% error 
c Rcst energy of line (keV) 
^Intrinsic width of line (eV) 

e Observed intensity of line (10~ 5 photons/cm 2 /s) with 68% error 

^Equivalent width of line (eV). 90% upper limit is shown if the K/3 line is very weak. 

9 Observed intensity of line (10~ 6 photons/cm 2 /s) with 68% error. 90% upper limit is shown if the K/3 line is 
very weak. 

h The increase of x 2 after removing K/3 line from the best-fit model 

l The increase of x 2 after removing Ka and K/3 lines from the best-fit model 

J The rest energy of line is fixed at the value obtained by the co-add spectra, since the K/3 line is very weak 
in this observation 



Table 4: The derived parameters from the spectral fits of the co-added XIS03 and PIN data (with error corresponding 
to 90% confidence level except flux and intensity). 



Sequence 
number 


Continuum 3, 

rb E cut-off c Rd 


Flux b 


E b 


Ka 

a h I b EW b 


E b 


K/3 
I b 


EW b 


X 2 /dof b 


A Xl b 


Axi b 


1 


1 5Q+0-03 jr+85 7q +0.35 
i - oy -0.03 ,0 -15 u '' y -0.32 


1 r i +0.1 


°- oyu -0.007 


38+ 16 2 02 +al3 95+ 21 


'• U8 -0.05 


3 4+0-9 

J -^-0.9 


i9±i; 


423.8/396 


4 7 


405.5 



o 

I 



a The cosine of the disc inclination was fixed at 0.87 
fc The meaning is the same as that in Table 3 
c Energy cut-off (kcV) 

^Reflection parameter of cold matter; R = £1/2tt 
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